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We study charge screening in a system of two dimensional nonabelian vortices, at 
_L \ finite temperature. Such vortices are generated after an SO (3) global symmetry group 

is spontaneously broken to a discrete subgroup Q 8 , where the latter is isomorphic to 

O ■ 

the nonabelian group of quaternions. Poisson-Boltzmann like equations are derived 

• i-H , 

for the various inter-vortex potentials, and the solutions to lowest order in a small 
fugacity expansion, are shown to behave much like those in abelian or Coulomb charge 
systems. The consequences for the phase structure of the system are discussed, where 
the fugacities associated to the thermal production of the various types of nonabelian 
vortices, are shown to play a key role. 
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1 Introduction 

Vortex defects in two and three spatial dimensions, have many interesting properties 
which have found a wide variety of applications in two notable areas, condensed matter 
physics [1], and cosmology [2]. Indeed it is well known [3] that the occurrence of vortices 
in 3-dimensional condensed systems such as superfluid helium IV, play an important 
role in understanding the physical properties of these systems. 

Vortex defects in any dimension are characterized by 7Ti(M), the fundamental group 
of the vacuum manifold M of the theory. The type of vortex defects that occur in for 
example helium IV are 'abelian' in the sense that they are characterized by an abelian 
fundamental group, isomorphic to the group of integers Z. In 2-spatial dimensions, 
abelian vortices have particularly simple interactions which allows one to write down 
the grand canonical partition function for thermal pair creation to all orders in the 
vortex fugacity in terms of the so called Coulomb gas model [4] . The pioneering work 
of Berezinsky [5], and Kosterlitz and Thouless [6], into their statistical mechanical 
properties, showed that a gas of such vortices underwent a novel kind of phase tran- 
sition at some critical temperature T c . A simple physical picture emerged where for 
T < T c , vortices and antivortices form a medium of bound pairs which subsequently 
dissociates into free vortices and antivortices for T > T c . These results when applied 
to approximately 2-dimensional systems such as helium IV thin films, lead Kosterlitz 
and Nelson [7] to predict a universal jump in the superfluid density at T c , which was 
later experimentally verified [8]. 

In a previous publication [9] , we presented the results of some preliminary investiga- 
tions into the nature of interactions of nonabelian vortices in a particular 2-dimensional 
model. In this paper we shall focus on the problem of screening of the forces between 
nonabelian vortices, again in the context of the simple model considered in [9] . Ordi- 
nary abelian charge screening is the mechanism behind the Kosterlitz-Thouless phase 
transition [1]. Although there are a number of ways of exhibiting this screening, for 
example by exploiting the map between the vortex gas partition function and that of 



sine-Gordon theory in 2-dimensions [10], perhaps the most illuminating is to derive 
a so called Poisson-Boltzmann (P-B) equation [11], satisfied by the linearly screened 
potential between a test vortex and antivortex placed in the gas. One of the goals of 
the present paper is to derive a P-B equation for the various inter-vortex potentials for 
the case of nonabelian vortices. In the model under consideration, there are basically 
2 species or types of vortices in the nonabelian system at finite temperature. We shall 
show that to first order in a small fugacity expansion, the P-B equations effectively 
become those of a pair of 'abelian' P-B equations , one for each vortex type. However 
because of the nonabelian nature of the vortices in question, the system is not simply 
reducible into two non-interacting 'Coulomb' systems to this order. In fact, because of 
the nonabelian fundamental group 7Ti(M) , we find that the coupling constants of the 
two Coulomb systems are constrained. 

The structure of the paper is as follows. In Section 2 we give a brief review of the 
2-dimensional field theory model [12] which describes nonabelian vortices associated 
with 7Ti(M) = Q 8 , and which was studied further in [9]. Vortices corresponding to this 
fundamental group may have some relevance to certain liquid crystals [13]. In section 
3, as an introduction to Poisson-Boltzmann equations and to fix notation, we show 
how to derive the P-B equation for a Coulomb gas, which is known to be physically 
equivalent to a system of abelian vortices. Sections 4 and 5 are concerned with deriving 
the same for the nonabelian vortex system defined in section2, and some conclusions 
are drawn about the possible phase structure of the system. 

2 A model of two dimensional nonabelian vortices 

Perhaps the simplest model [12] in which nonabelian vortex defects occur, is that in 
which spontaneous symmetry breaking occurs in two spatial dimensions, with an order 
parameter <3> characterizing a system whose total energy E we may choose to be 



E = f d 2 x - Tr [g ab d a $d b <S>] + V($) (2.1) 



V($) = -Tr$ 4 + -(Tr$ 2 ) 2 

+ ^Tr$ 3 -^ 2 Tr$ 2 (2.2) 

where $ is a scalar field transforming in the five dimensional representation of the 
symmetry group G = SO (3), i.e. it is a traceless 3x3 symmetric matrix. 

When p ^ 0, V has three isolated minima and the unbroken symmetry group is 
U{1). This case is not particularly interesting since there are no stable vortex defects 
produced. This is a result of a topological theorem which states that if the first ho- 
motopy group 7Ti(G/H) is trivial, where H is the unbroken symmetry group, then the 
vortices are unstable, and indeed in this case tti(SO(3)/U(1)) ~ 1. If p = however, 
50(3) is broken down to a discrete subgroup K isomorphic to Z 4 the additive group 
of order 4, and in this case there is a degenerate set of minima of V which lie on the 
ellipse tp\ + (fl + (pi<p2 = / i o 2 /(-^ + A ) , where the latter two fields are the independent 
eigenvalues of the diagonalized field $. The fundamental group of the manifold of 
vacuum states of the theory defined in eqs.( [2.1|) and ( |2.2|) , is given by 

7nOSO(3)/IK)~Q 8 (2.3) 

Q 8 being a nonabelian discrete group of order 8, (isomorphic to the group of quater- 
nions). It is generated by elements i, j, k, —i, —j, —k and -1, with 

f = f = e = _ t 

ij = k , jk = i , ki = j (2.4) 

where i, j, k define a basis of quaternions. Hence there are three types of nonabelian 
vortices in this model corresponding to i,j and k, ( the elements —i, —j and —k refer to 
the corresponding anti- vortices, whilst —1 actually defines a Z2 vortex which is abelian 
in nature since —1 commutes with all the other elements in Q 8 ). The vortices which 
are produced by this symmetry breaking are guaranteed to be stable by topological 
arguments. 

In ref.[9], the following ansatz was used to minimize (numerically) the energy E 

$(r,0) = £(fl) W)^), (2-5) 



where the diagonal matrix $dia g = ( { Pi,f2,—(^Pi + ¥2))- The SO (3) group elements 
Q{9) satisfy the nonabelian vortex boundary conditions 

g(9 + 2n) = Q{9)h (2.6) 

where h is an element of the group IK, and (r, 9) are polar coordinates centered on the 
vortex core. Explicit solutions for the functions <pi,<fi2 and Q pertinent to describing 
vortices of various types, were given in [9]. Moreover, it was also shown that the ansatz 
in eq.(j2.5|) was stable to perturbations that preserve the boundary conditions, eq.( [2.6|) . 



In this sense, lowest energy vortices are correctly described by eq. fl2.5|) . 

For the energy of the vortex to be finite, $(r, 9) must tend asymptotically as r — > 00, 
to a minima of the potential V(§). In addition, $(r, 9) should vanish sufficiently 
quickly as r — ► 0. Classical field configurations $ with these properties were obtained 
by numerical minimization of the energy E and the reader is referred to [9] for further 
details. 

Having mentioned the stability of vortices described by eq.(|2.5p to perturbations 
that maintain the homotopy classes, one still has to be cautious before concluding that 
such vortices will be stable if produced thermally in a system at finite temperature. 
This is because dissociation can occur between one vortex type and another, depending 
on the relative size of the chemical potentials involved. From a group theoretic point 
of view, this process of dissociation can be understood simply as a consequence of 
group multiplication, e.g. the relation ij = k in eq. fl2.4p implies that a k type vortices 
could decay into an i and j type. In [9] it was shown that the chemical potentials 
controlling the concentration of the various types of vortices, depends on the single 

real parameter a = , which essentially fixes a point on the ellipse of vacuum 

states in this model. As the value of this parameter is varied, either one of i, j or k 
type vortices are unstable to dissociation into the other two. For definiteness, we will 
restrict a to be in a range such that k types are unstable to dissociation into those of 
i and j. (In passing, it was proved in [9] that the abelian vortices corresponding to 
the element —1 always dissociates into vortex-antivortex pairs for any choice of a, so 



we will ignore these in the remain of the paper). To summarize, kinematic constraints 
simplify the study of nonabelian vortices at finite temperature, at least in the model 
under consideration, by allowing us to focus on just two types of vortices corresponding 
to particular non-commuting elements of Q 8 . 

It is the purpose of this paper to try and understand the screening properties of the 
(i and j ) type nonabelian vortices, by deriving Poisson-Boltzmann (P-B) like equations 
for the corresponding inter- vortex potentials. Since the P-B equations are more familiar 
in the context of screening in systems of abelian vortices, in the next section we shall 
review their derivation, and the consequences they have for Kosterlitz-Thouless type 
phase transitions. 

3 Coulomb gas system 

It is well known that a system of abelian vortices in 2-dimensions can be mapped 
onto that of the so called Coulomb gas [4]. Abelian means that the fundamental 
group associated with spontaneous symmetry breaking is itself abelian. Such vortices 
occur when U(l) symmetry is broken e.g. in superfluids, superconductors etc. [1]. 
Because of the mapping between systems of Coulomb charges and abelian vortices, the 
notions of charge and vorticity are interchangeable, and we shall often use both sets of 
terminology in what follows. Consider then, a Coulomb system of N charged particles, 
with charges ±q, and let N + ( AL ) be the total number of positive (negative ) charges 
with N = N + + AL. Moreover, let us restrict ourselves to neutral Coulomb gases, i.e. 
those in which N + = N_. The total energy of such a configuration is given by 

N 1 

h n = £ ~Qa<lbU{r ab ) -Nfi (3.1) 

where in eq.( |3.1|) , U(r a b) is the electrostatic energy of charge q a due to charge q b and 
fab — | r o — r&| is the distance between those charges. The chemical potential \x is a 
sum of the electrostatic self-energy U(0) and non-electrostatic contribution, E c or the 



core energy. 

p= -.[lu{0) + E c ] (3.2) 

The partition function of the system is 

oo -I AT . t2 

^=e o m e n/^-"* (3.3) 

where 7r£ 2 is a phase-space division factor and £ is a length scale typically the size 



of a Coulomb particle. In eq. (|3.3| ) and subsequently, we set the Boltzmann constant 
ks = 1. For a system of interacting abelian vortices, £ corresponds to the core 
radius. The restriction J2Qa — on the sum over charges in eq.( |3.3| ) follows from the 
assumption that we are considering a neutral system. Introducing two infinitesimal 
test charges 5q + , Sq- into the system at fixed positions r + and r_, we can approximate 
the partition function of this modified system as follows 

Z(5g,r+,r_) = 1+ fd 2 re^+^ (r)/ V(r-(r + -r_)) 

J 7f£ 2 
r dj^ [l6q + q_V L (T-r + )+±V L (0)+E c ]/T-[±Sq-q-.V L (r-r-)+±V L (0)+Ec]/T 

J *e 

+ ... (3.4) 

where in eq.( |3.4j ), we have only kept terms at most linear in the fugacity z = e~ Ec ^ T , 
which is taken to be small. 

Since the two test charges introduced into the Coulomb system are not thermally 
created, there will be no chemical potentials appearing in the second term of eq.( |3.4| ). 
Vl(t) is the linearly screened potential at r which in the absence of background 
charges, is simply the unscreened Coulomb potential depending logarithmically on 
distance between charges. The interaction between the test charges Sq + and Sq- and 
the background charges q + and g_ which are thermally created, depends on the core 
energy E c and the self-energy Vl(0). We define 5q + = —5q_ = —Sq, q + = — g_ = q 
and the density of positive (negative) background Coulomb gas charges n F as 

«?' = -—7, = -^-' i,0,/2r (3-5) 



then we may rewrite the partition function as 

Z(6q, r+, r_) = 1 + f d 2 v e ^ Sq)2yL{r)/T 5 2 {v - (r+ - r_)) 



d 2 r 



rl 2 Y- 

11 l -y q qV L (r-r+)/T +±6qqV L (r-r-)/T 



+ n 

" J 7T4 Z 

D + UqqV L (r-T+)/T -UqqV L (r-T-)/T 



n £2 



(3.6) 



To obtain the linear screening potential at the point ro = (r + — r_), we apply techniques 
which are familiar in field theory. First introduce a source term or test charge into the 
generating function or partition function (as we have done above). Then, differentiating 
the new partition function with respect to these sources gives us a correlation function 
which is proportional to Vl, 

5qV L (r + -r ) = T dXnZ ^ r +' r ~) - T dZ(5g,r + ,V-) 

d 5q Z(5q, r + ,r_) d5q 

Thus, 

5qV L (r + -r_) = || 5q V L (r + - r_) e |(M 2 ^(r + -r_)/T 

+ IJlZ fd 2 r [-V L (r - r+) + V L (r - r _)] e -|*ff«V £ (r-r + )/T + §* M V i (r-r_)/T 
+ iHZ. /> r [V L {T - r + ) - V L (r - r _)] e ^qqV L (r~r + )/T -l5qqV L (r-r„)/T + \ 

(3.8) 

Since we assume small fugacity z and shall ignore nonlinear terms in Sq (by taking 
Sq sufficiently small), one can approximate 1/Z by unity (Z = 1 + 0(z)). To obtain 
the Poisson-Boltzmann equation [1, 11], we apply the two dimensional Laplacian V 2 
with respect to the relative position of the test charges r = (r + — r_) to both sides of 
eq.d; 



8qV 2 T0 V L (y + - r_) « 5q V r 2 V (r + - r _) e ^) 2 ^ + —)/ r 

+ qnFJ d 2 r v2J-y ( r - r + ) + V (r - r_)] e -^«« v *< r - r +>/ T +±5 qq V L (r-r_)/T 



+ ^E. Jd 2 T V 2 JV (r - T+) - V (r - r _)] e |*9«Vb(r-r + )/T -l Sqq V L (r-r.)/T + 

(3.9) 

where Vo(r) is the lowest order (unscreened) approximation of V^(r). Since Vo(r) is the 
interaction potential of charge particles in two spatial dimensions, it is proportional to 
the logarithmic function, so that 

V 2 T0 V (r ~ r+) = -2na6 2 (r - r+) (3.10) 

where a is a constant with dimensions of temperature, and is the effective coupling 
constant of the Coulomb gas. Thus we obtain the Poisson-Boltzmann equation as 



SqV 2 V L {r + - r_) = -27ra5q 5 2 (r + - r_ 

2naqn F 1 

e -^dqq[V L (r + -r-)-V L (i))\/T {3-H) 



2naqn F ^ Sqq [V L (r + -r^)-v L (o)]/T 

27TO q Tip ^l Sqg [ VL (r + -r-)~V L {0)}/T 



2 

where n F = n~^ + n F is the number density of free charges. Next, we redefine 
Vl(t + — r_) as the separation energy between a test charge and background charge, 

E Sep = V L (r+ - r_) - V L (0) -> Vi(r + - r_) 

where Vl(0) is the energy of zero separation between test charge and background 
charge. Finally, dropping the ~ (tilde) and expanding the exponentials to lowest order 
in 5q we obtain the Poisson-Boltzmann equation corresponding to an abelian (in the 
sense defined earlier) Coulomb gas 

r 2 

8qV 2 V L {v + -rJ) = -2na5q5 2 (r + -r_)+ — ^ nF V L {r + -T-) (3.12) 



Solutions to eq. (|3.12j ) have been discussed in the literature [1,11] in the context of 



charge unbinding in the Coulomb system. It is clear that when n F ^ 0, the potential 
V F is screened, and decreases exponentially with the separation between a pair of op- 
posite charges (or vortex-antivortex pair in the case of abelian vortices). Consequently, 
the pairs become unstable to dissociation into free charges, and so 'unbinding' of the 

8 



charge or vorticity occurs. This in turn implies of course, a non-vanishing density of 
free charges rip so is a consistent physical picture. The alternative possibility is that 
rip = 0, in which case the P-B equation for Vp implies that the latter is logarithmically 
decreasing with the pair separation, giving rise to a medium of bound neutral pairs. 
The question of whether rip is vanishing or not depends on the temperature T of the 
system, i.e. rip = if T < T c , and rip ^ if T > T c , where the critical temperature 
T c = — (1 — 2nz) define a line of values as z is varied [5,6] . T c marks the position of 

the so called Kosterlitz-Thouless phase transition of the Coulomb gas [6] . Perhaps the 

a 
simplest way of obtaining the leading order ( in z ) approximation of T c , (i.e. T c = — ) 

is by the method of self-consistent screening [7] . One introduces the so called screening 

length uj(z, T), which for an infinite system is related to the density of free charges rip 

by 

T 

n F (z,T) = (3.13) 

7iauj z {z, 1 ) 

where it is clear from eq. (|3.13|) that uo depends on the fugacity z and T. At the same 



time, since the linearly screened potential at the origin Vl(0) = hi(y), eqs.(|3~5|) and 



? 



(|3.13|) together provide a consistent set of equations determining the value of uo . There 



are two self-consistent solutions namely for T > T c , oj is finite (and Vp exponentially 

decays ) whilst for T < T c the screening length is infinite and the potential remains 

a 
logarithmic, with T c = — [7]. Alternatively, one can deduce the dependence of A on 

z and T, by utilising the well known map of the Coulomb gas onto the sine-Gordon 

model [10], and exploiting the known renormalization properties of the latter. This 

method also gives the next to leading order form of T c given above. In section 5, we 

shall return to the problem of solving the Poisson-Boltzmann equation for systems of 

nonabelian vortices or 'charges', which we shall introduce in the next section. 

4 System of nonabelian vortices. 

We now wish to turn our attentions to a system of two dimensional nonabelian vortices 
of type i and j, (i.e. field configurations satisfying the ansatz of eq.(|2.5|) with twisted 



boundary conditions as in eq.( |2.6| ) with h = i and j respectively.) To construct the 
Poisson-Boltzmann equation for the nonabelian vortex system, we first have to define 
the infinitesimal charge of a nonabelian vortex. From our ansatz eq.( [2.5|) , the field 
configuration $ depends on the SO(3) group elements Q a , which represent rotations 
about the x, y, or z axes. Therefore, the infinitesimal charges are defined correspond- 
ingly as rotations through an infinitesimal angle e9 about those axes. For example, 
the classical field configuration $ of an infinitesimal test charge of type i is 

$ei = Gei®diagG7? (4.1) 

where 



(I \ 

cos(^) sin(^) 
VO -sin(^) cos(^)7 



(4.2) 



with similar definitions for the vortices of type j and k. One can compute the infinites- 
imal nonabelian vorticity or nonabelian charge of these classical field configurations, 
by for example calculating the form of the self energy in each case and comparing it to 
the abelian case. One finds 



5qi = 


-8q-i 


= -7|(^i + 2 ^2) 


5 qj = 


-8q-j 


= 7^(2^1 + ^2) 


Sqk = 


-Sq- k 


= ^1-^2) 



(4.3) 



with infinitesimal parameter e. The corresponding finite nonabelian charges are 

1 



q-i 


= -q-i 


= -y|(^l+2<^2) 


Qj 


= ~ <i-j 


1 . 
= -^(2^1 + ^2) 


qk 


= -q-k 


= ^1-^2) 



(4.4) 

It should be noted that the "charges" defined in eq.( |4.4| ) have dimensions of ^/energy, 
compared to the dimensionless ones introduced in section 3. 

10 



Having discussed how to define the notion of infinitesimal charge for nonabelian 
vortices, we are now in a position to consider the 'nonabelian' charge-unbinding scheme 
via the Poisson-Boltzmann description. Let (iV_j , JV_j) iV+j , N + j be the total number 
of % and j type (anti) vortices present. As usual we define JVj = N + i + N^i and 
Nj = N + j + N_j and again restrict ourselves to neutral systems i.e. to configurations 
having equal numbers of vortices and antivortices of the same type, so N +i = N_i 
similarly for j types. 

Under these circumstances, the partition function of the system is 



oo oo i Ni Nj 

z = e E^nn E E 

Ni =oN j= o iV - J V ".=i« 3 =i goi ;E s . 9oi =o ga ,;E a . 9o .=o 
x /%/%e-^^)/ T (4.5) 

Generally, one can write down the hamiltonian of the nonabelian system as 

Ni l Nj l 

H(N h Nj) = Y, 2 qa ' q ^ Uii ^ ~ n ^ + E 2 qajqbjUjj ( Taj ~ n ^ 

N u Nj 1 

+ E 2 qaiQbjUij( ' rai ~ n ^~ ^ iNi ~ ^ jNj ( 4 - 6 ) 

ai, bj=l 

The quantity Ua(r a . — r b .) is the interaction energy due to the charges q ai , q bi of the i 
type vortices amongst themselves, with similar definitions for Ujj and Uij. The chemical 
potentials fii and fij are 

H = -[\u n (0) + Ei] (4.7) 

where E l c and E 3 C are the corresponding core energies of % and j type vortices respec- 
tively, explicit expression for which may be found in [9]. We define the fugacities 
corresponding to the thermal creation of vortices of type % and j respectively 

Zi = e~ E ^ T , Zj = e- Ei ' T . (4.8) 

Now we introduce two infinitesimal test charges bqi and Sq-i into our system at the 
points r ei and r_ ei where we remind the reader that only the vortices of type % and j 
need be considered if we restrict the parameter a in the range (—2 < a < — |) ). 

11 



Then, expanding the partition function denoted by Z^Sqi, r ei , r_ ti ) for the modified 
system in the presence of such test charges to linear order in the fugacities (which are 
again taken to be small), we obtain 

Zi{S qi , r ei , r_ eJ ) =1+Jd 2 r g-^-^M/T^ _ ^ _ r _ ej)) 

+ Zi [ (fr e -lEHi(r-r ei )+E- Hi (r-r- ti )]/T 

+ Zi f d 2 r e~ lE "- l{r ~ r " )+E -"- tir ~ r -" )]/T 
+ z- [ d 2 r e -l E *ii( r - r *i)+ E -tij( r - r -*i)]/ T 

+ z- [ d 2 r e~ [E "- j{r ~ r " )+E -"- j(T ~ r -" )]/T 

+ (4.9) 

where the subscript i on Z indicates the partition function in the presence of an in- 
finitesimal i type vortex and V" is the linearly screened potential between a vortex- 
antivortex pair of type i. The interaction energy E e u(r — r e ») = SqiqiUu(r — r 6 j) with 
similar definitions for the other interaction energies appearing in eqs. Q4.9p in terms of 
the interaction potentials Ujj and Uij. According to our ansatz for the classical field 
configuration $ of a nonabelian vortex configuration, is 

Q = G$d iat ,G- 1 (4.10) 

and the energy of such a configuration can be approximated by the kinetic energy term 
in eq. (pTT|) 

E « I d 2 r Tr 



g ab d a $d b <f>\ (4.11) 

Therefore, it is clear that the unscreened interaction energy E® u (ti — r ti ) between a 
background charge qi and infinitesimal charge 5 qi at points r, and r e j is given by 

/ b d a <S> eU d b <S> eli 
J d 2 v^2e 2 {ip x + 2 if 2 ) 2 [W r A{v - v el )] 2 

+4e(<^i + 2<p 2 ) 2 V r A(r - r d ) ■ V r A(r - n) 

12 



E&fa - r et = / dh Tr 



+2(^ 1 + 2^ 2 ) 2 [V r A(r-r i )] 2 } 
= 2(6q i ) 2 VJ i (0) + A5q t qtVfin - r el ) + 2( ft ) V o "(0) (4.12) 

where 5qi — eqi, A(r) = Im [ln(re j61 )] and 

V^iri-Tei) = JdhiVrAir-r^-iVrAir-r,)} 

(4.13) 

is the unscreened interaction potential between two i-type vortices. Furthermore, 

Kf(0) = /cfrpvVKr-r,)] 2 

= Jd 2 r[V r A(r-r ei )} 2 (4.14) 

are the self-energy of the background and test vortices at the points r and r e j re- 
spectively. Since the test charge <% is not thermally created, we should subtract it's 
self-energy term 2(6qi) 2 Vq u (0) from the right hand side of eq. (|4.12j ) Henceforth, we 



shall take for the form of the interaction energy including the effects of screening, En, 
the following 

E it = Adq^Vl 1 ^ - r €l ) + 2q 2 V?(0) (4.15) 

( This is analogous to the replacement of the unscreened logarithmic Coulomb potential 
by Vl in the previous section.) Similarly, 

E. ai = -46q t q i Vl t (r l -r. el )+2q 2 Vl l (0) 
E €l _ t = -A6q l q i Vl i (r i -r €l ) + 2q 2 Vl i (0) 
£_«_, = 45q iqi V?(r t -r_ el )+2q 2 V?(0) (4.16) 

Next we shall consider interaction energies between an i type test vortex and a 
background vortex of type j. This now involves 2 possible orientations of the group el- 
ements Q e i and Qj, since the latter elements are non-commuting, but a straight forward 
calculation gives 

^tij "' yeiyj^diagyj =^ei ' jti =J j =4 ei* diag =J ti ^j 

13 



E eij 



-{ / d 2 v Tr g ab d a $ elJ d b $ 



3 



+ / d 2 v Tr 



g ab d a $ jei d b $ j€i } 



J d 2 r{ ^2 e 2 (<^ + 2 ip 2 ) 2 - -eV 2 (2^i + <^ 2 ) sin 2 [A(r - r,)]) [V r A(r - r e 
2(2 v?i + ^ 2 ) 2 - 2 e V 2 (2^i + ^2) sin 2 [A(r - r ei )]\ [V r A{r - r,)] 2 } 



(4.17) 



whilst for a background anti- vortex of type j we have 



E ei~j 



GeiGj ®diagGjG ei} §-j ei =Gj Gei^diag G ei Gj 



d 2 r Tr 



g ab d a <f> eij d b $ eij + / d 2 r Tr g ab d a ^ jei d b ^j 



,06; 



£, 



«j 



(4.18) 



From eqs. (|4.17f ), (|4.18|) one can see that in fact, there are no interaction terms between 
the two different type of charges Sqi and qj, which would have corresponded to terms 
linear in e. The terms present can be interpreted as the self energies of the % type 
test charge and the j type background charge. It should be noticed here that the self 
energy of the test charge Sq^ has been modified by the sin 2 terms. Since the test charge 
is not thermally created, all such self energy terms should be subtracted from E € ij as 
stated previously. 

As is shown in appendix A, the modification terms themselves do indeed take the 
form of divergent self energies of the test charge Sqi , so we are justified in dropping 
them from the expressions for the interaction energies. Hence finally we have 



E, 



IIJ ^2q 2 Vi J (0) = E- eij 



E, 



€l-J 



E_ 



ei-j 



(4.19) 



Similarly, one can check that 



E 



E_ 



ejj 



<:.!.) 



48q jqj Vi J (r-r ej )+2q 2 Vi J (0) 
-U mj Vl J (r - r. ej ) + 2q 2 V[ j (0) 



E 



E_ 



c-J-3 



Q-3 



-48q jqj Vi J (r 
4Sq iqi Vl J (r- 



- r ej ) + 2q 2 V[\0) 
v. ej ) + 2q 2 Vl\0) 



(4.20) 
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and 

E eji = 2(gVi i {0) = E_ eji = E ej _i = E_ eJ _ t (4.21) 

The density of positive(negative) background nonabelian charges is defined as follows 

(+) e (-4V«(0)-.E*)/T 



'i 


*e 




ne 


(+) 


e (-wl J (o)-Ei)/T 


J 


ire 




- z i p -4V^'(0)/T 



(4.22) 



ire 

Substituting the expressions for the total energies E of the various configuration as 
calculated above in terms of the linearly screened potentials, we obtain for the partition 
function of eq. ( fX~9[ ) 



Zi(S qi , r eu r„ el ) = 1 + J d 2 v e -¥^->v?(r)5\r - (r el - r_«)) 

+ z~ [ d 2 r e -[- 4 ^<?^2 l ( r - r «)+ 4 %<?^2 l ( r - r -«)+ 4 9?v«(o)]/r 

+ zfJd 2 re^ v £ 3 ^ T 

+ zj [d 2 re^^ J ^ T + ... (4.23) 



Rewriting the partition function in terms of the free charge densities, we obtain 
Zk = 1 + j d 2 ve^) 2 W5 2 {v-{v el -v^)) 

+ Tlf f d 2 T e -4-5^g4V2 l (r-r El )-^(r-r_„)]/T 

+ n~ J d 2 V e -^mhV-(r-r a)+ V-(r-r.^/T + ^^ 

5 P-B equations for nonabelian vortices 

Before we compute the Poisson-Boltzmann equations of the system, we shall consider 
the relationship between the two different types of nonabelian vortex charge. The 
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u 



-1 0^ 
0-10 
1 
/-l 
1 

V 0-1, 

1 0^ 

1 
1, 



Qi = ug.u- 1 



9i = ug.u- 1 



G-i 



Qj 



Qi 



U < <5>diaq > U~ 



¥>1 





/<Pl 



Vo 






'■p2 




<r2 



V>i 








-(cpi +<p%) t 



-((ft + (fi 2 ), 



-(^1 + ^2)- 



Table 1: 

energy density (again dropping the potential term as explained earlier) of a single i 
type vortex is written as 



Si = Tr [g ab d a ^d b % 

which, up to core corrections, can be approximated by 

1 



Si 



2r : 



:((<p 1 + 2<p 2 )V r A[r-r 1 ]Y 



(5.1) 



(5.2) 



at distance scales greater than the typical vortex core size £. In eq.( |5.2| ) we have 
considered a single vortex of type i located at the point ri and the field configuration 
$i is given by 

^i t'i^aiagt'— i 

as described earlier. 

Similarly, for a single j-type vortex at the point r 2 , 



S 3 = Tr 



g^dafydb® 



j^b-^j 



Sj = —^((2^ + (p 2 )V r A[r - r 2 ])2 



2r 



(5.3) 
(5.4) 



16 



From table 1, it is clear that we can transform Qi to Qj by using the constant transfor- 
mation matrix U, so we can rewrite $j in terms of Qi as follows 

$. = UQ-U-^^U-^Q^U (5.5) 

Thus we could have obtained the energy density of the j-type vortex, by substituting 
using eq.( |5.5| ) into the self-energy of an i type vortex. Since we know that U is unitary, 
by the property of the trace which is invariant under the similarity transformation 

Tr [B] = Tr [UBIJ- 1 ] 

we then conclude that 

Tr 



g^d&dbQj = Tr g^d^d^ (5.6) 

which is evident from comparing Si with Sj . It follows that the infinitesimal charges 
Sqi and Sqj transform into one another by simply interchanging tp% and tp 2 in $diag- 
This is sufficient to show that both infinitesimal charges Sqi and 5qj depend on the 
same infinitesimal rotation parameter e. 

To find the linear screening potential of the system with both types of infinitesimal 
charge 5qi and 5qj present, we simply use the chain rule by differentiating the partition 
function first with respect to e and apply the following relations 

Sqi = -j=((fi + 2 ip 2 ) 

dSqi 1 

—=(Lp 1 + 2Lp 2 ) = q t 



Se y/2 

s qj = j 

d 5qj 1 



H = -^(2^! + ^ 2 ) 



(2^ + ^ 2 ) = q< (5.7) 



Se ^2 

For later use, we shall need to consider the case when infinitesimal test charges of type i 
and j are simultaneously present in the system. For this system with two different test 
charges 5qi and Sqj at r ei and r t j, the lowest order partition function can be reduced 
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to 



;%(%,<%, r d ,r 6i ) = l + Jd 2 re-^ 5 v v ?U/ T 8\v-(r ei -r ej )) 

+Zi f d 2 r e -{E e u+E eji ]/T 

+Zi f d 2 V e -[Eei-i+E ej -i]/T 

+ Zj f d 2 r e-l^ij+E^/T 

+ Zj fAe-^ +E -'/ T + ... (5. 



where the subscript ij on Z in eq. (|5.8|) denotes the fact that both i and j type test 
charges are present. Substituting the expressions for the total energy E from the 
previous section, we get 

Z ii (^ i ,%,r d ,r 6J .) = 1 + J d 2 r e^ 5 ^v^ r) /T 5 2 {Y _ {Yei _ rei)) 

+ nf /> re -«**V*(r-r ri )/T 

+ nrJd 2 re 4S ^ v ^- r ^ T 
+ n+ j ' d 2 ve i&q ^ v " (r - r ^/ T 

+ nj J d 2 V e -4%%Vf(r-r ej )/T + __ (5g) 

The linearly screened potential at (r ei — r eJ ) is given by 

W(r« - r ei ) = T 9 ^- = 1-™*L- (5.10) 

and 

Now we are in a position to construct a Poisson-Boltzmann equation for the linearly 
screened potentials between various test charges. In the first case, two test charge will 
be taken to be of the same type, whilst in the second case we will take different types. 
For the same type of test charges, say 5qi and 5q_i, the linearly screened potential 
at (r ei - r_ £i ) is given by 

^• ( ,, .,,„). r|^=l^ (5,2) 
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Since we use the two infinitesimal charges Sqi, Sqj, we can ignore the higher order terms 
and therefore approximate 1/Zj by 1 (Zi = 1 + 0(zi, Zj) + 0(5qf) + 0(6qf)). Hence 
we find 



lrx„.\2 



5 qi Vl\v el -v_ et ) = 5 qi Vl\v el -v_ ei )e^ d ^ ^ 



(%) 2 V»(r ri -r_ ri )/T 



4q t n\ 



T 

Aqi n~ 



fd 2 T [V?{r - Tei) - Vf (r - P_ei) ] e" 4 **® [V^*(r-r Bi )-y«(r-r_« i )]/T 
fd 2 r [-Vt\v - r d ) + V[\r - r _ ei )] e -*S qiqi [-V?(r-r ei )+V?(r-r_ ei )]/T 

(5.13) 

To lowest order V[ l (r) = Vq 1 {t) + 0(zi, Zj) where the unscreened potential V^*(r) varies 
logarithmically with distance. Applying the Laplacian operator V 2 to both sides of 
eq. (|5.13f) where r is the relative position of test charges r = r €i — r_ ei , we get 



5 qi Wl o Vl\r el -r_ el ) = Sq^l V™(r ei - r_ d ) e*<**> 3 ^(^—^ 



^ nj 



T 



fd 2 rV 2 ro [Vt\v - r ei ) - V?(r - r _ e .)] e -4«v»(r-r ei )-v«(r- r -ei)]/T 



4q t nl 



/d 2 rV 2 J-Vf (r - r €i ) + V?(r - r _ ei )] e -^ qi [-vi\r-^)+vi 



T 

(5.14) 

Since, as mentioned, the unscreened potential V^ M (r e j — r__ e j) is proportional to the 
logarithm function, it follows that 

V r 2 V eii (r - r d ) = -2vr5 2 (r - r d ) (5.15) 



Hence, using eqs.( [5.14! ), fl5.15| ) we find a Poisson-Boltzmann equation satisfied by 

^V r 2 ^(r e , - r_ e ,) = -27r^5 2 (r eJ -r_ eJ ) 

- Svrg, nT e -««<«(r. i -r_. i )-v*(o)]/T (5 . 16) 

Again as in the abelian case we redefine 

Vr(r ei - r_ ei ) - V?(0) - V*('ei - r_ei) (5.17) 
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We get finally, (dropping ~ (tilde) ) 

V 2 Vf (r« - r_ ei ) = -27T 5 2 {y a - r_ ei ) - ^ ^' H V?{v ti - r_ ei ) (5.18) 



327r a,; n; . 



where nf + n~ = n« and ctj = g^ 2 is the analogue of the quantity a introduced in 
section 2. 

Similarly, for the two test charges of j and —j type, we can derive a P-B equation 
for the corresponding linearly screened potential V[ J (r e j — r_ eJ ) 



V r 2 ^'(r y - r_ e ,) = -2vr 5\r eJ - r. ej ) - °—^Vl\v eJ - r_ ej ) (5.19) 



with n^ + rij = rij and a.,- = g^ 2 It should be remembered that in deriving the P-B 
eqs. (|5.16|) and ( |5.19| ), we have used the property, that V^ is vanishing at lowest order, 



as discussed in section 4 immediately following eqs.( 4.18|) . This is why V^ plays no 



role in the P-B equations derived above. A consistency condition on this property of 
the screened potential can be derived by considering a P-B equation for V^ itself, and 
showing that the solution to lowest order is indeed vanishing. To derive such a P-B 
equation, we consider a simultaneous configuration of different types of infinitesimal 
test charges Sqi, Sqj. From ( |5.10|) we obtain 



6qiV L J {r ei -r ej ) - 



q% de 

- 5 qi V[ j {r ei - r ej ) e-3***«* v^-r^/r 

- ^}L />r Vf(r - rei ) e ~^mViHr~r„)/T 

qj J 



Qj J 
+ Ml£ [d 2 r V[Hv - r„.) e -^*^(*-^>/ T 



Qj 

,2, 



A Qj n o 



d 2 r yjj^ _ v ^ e A5 mi Vi\v-v ei )/T 



To lowest order in Z{ and Zj, 



Qj 
+ ... (5.20) 



Vl\r) = V*\r) + 0{z u z 3 ) (5.21) 
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where as we have shown in Appendix A, Vq (r) is in fact a divergent self energy con- 
tribution to the test charge, independent of r. Applying the Laplacian operator Vp to 



both sides of eq. (|5.20|) , where r is the relative position of the test charges r = r ei — r e j, 
we obtain 



8 qi Vi V?(r ei -r ej ) 



5q t Vly L \r„ - r„) e^** 8 * ^(r E1 -r ej )/T 
Aq 2 nf 



+ 
+ 



Aq 2 nj 


qj T 

Aq 2 n+ 


qj T 
Aq)n~ 



qj T 



d 2 T 

d 2 v 
d 2 v 

dh 



VXr-r« 



V 2 ro V?(r-r e 



V^f(r-r e ,; 



VJoVjftr - r*J 



46q m V»(r-r ei )/T 



a 4Sq m Vl^r-r el )/T 



7 ASq jqj Vi J (T-v ej )/T 



46q jqj Vl>(r-r ej )/T 



+ ... 



(5.22) 



From the properties of the lowest order potentials V™ , Vq 3 and Vq 3 discussed above, 
the P-B equation becomes 



SqNi Vl 3 (r ei -r 



<:.!) 



wqt n i e -4% gi y«(o)/T 






v«(o)/t 



_ M( Jj H 3 c 48g iqi V?>(0)/T 



+ 



9j 

^% n 3 c -46q iqi V?H0)/T 
1j 



(5.23) 



For neutral system, we know that nf = n i , then expanding the exponential terms out 
to the lowest order in fugacities, we get the final form of the P-B equation for V^ 3 



V^(r el -r y ) 



32iraiqi n % u 



qj T 



VL l (0) 



32 ™3 n 3 V 33(ty 



(5.24) 



A similar equation can be derived for V[ l which is given by the r.h.s. of eq. (|5.24 ) but 
with the subscripts i and j interchanged. Notice that an important difference between 
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eq.( |5.24| ) and the previous two P-B equations, ( |5.18| ) and ( |5.19| ) is the absence of the 



(fugacity independent) delta function term on the r.h.s.. Indeed the r.h.s. of eq.( 5.24 ) 
is simply constant in the variable ro 3 . Before moving on to discuss the solutions and 
consequences of these equations, it is clear that the role of the potential V[ J can be 
ignored, to the order in which we work. 

A comparison of eqs. (|5.18|) and (|5.19| ) to that of the P-B equation derived in section 
2 for the abelian Coulomb gas, shows that the potentials V[ l and V[ J are essentially 
Coulombic in nature with the free charge density and fugacity np and z replaced by the 
same quantities with subscripts i, j appropriate for the system of nonabelian vortices. 
Consequently, the solutions of eqs. (|5.18| ) and ( |5.19| ) are exactly of the type described 



in section 3, with the replacement of fugacity z by Zi , Zj and screening length uj by 
Ui(zi,T) and Uj(zj,T). Furthermore we see that there will in principle be two lines of 
critical temperatures T % c and Tj! given by 

n> ■ 
Tl = ^(1-27^) 

Tl = -f(l-2ir Zj ) (5.25) 

It is perhaps worth commenting at this point, that the two critical temperatures given 
in eqs.( |5.25| ) for vanishing fugacities Z{ and Zj , are in agreement with those one could 
have obtained through heuristic free energy arguments first discussed in [6]. If we 
imagine increasing the temperature of the system from T = 0, eqs.( |5.25| ) indicate 
that either the % followed by the j type vortices dissociate or vice versa depending on 
the relative magnitude of the critical temperatures in eq.( |5.25 ). Note that since the 



quantities c^ and <x,- can be re-expressed in terms of the vacuum values of the fields (pi 

and <£>2 

1 2 

oti = -< (fi + 2ip 2 > 



3 Thus the vanishing of Vfj found previously to lowest order in fugacity is consistent with the P-B 
equation for V? 
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There is a constraint on ctj and a, since the vacuum values of (pi , tp 2 must lie on the 
ellipse ip\ + ip\ + y?iy?2 = f^o 2 /^ + A ) . This constraint in terms of an and a, is 

3„2 

on + aj - ^a~a~ = + ° (5.27) 

which shows that the values of ■s/oTi, JoTj also parameterize an ellipse. It is clear from 
this that there exists an ellipse of critical temperatures T* and T-l, as the fugacities z% 



and Zj are taken to zero, 



3/x 



2 



It + It - yjTiTi = ^f^ (5.28) 

This curious phase structure is clearly a consequence of the degeneracy present in the 
vacuum manifold which as we have seen in section 2, is an essential feature of the 
nonabelian vortices considered in this paper. 

In conclusion, we have investigated in this paper lowest order screening in a system 
of nonabelian vortices in 2-dimensions, through a detailed study of Poisson-Boltzmann 
like equations for the various inter-vortex potentials. Remarkably, we have seen that 
such equations are similar in form to those of an abelian system of vortices , the latter 
being physically equivalent to a Coulomb gas. This feature might be a consequence 
of the simple model of nonabelian vortices we studied and perhaps also to the fact 
that we considered only the lowest non-trivial order in a small fugacity expansion. 
It would certainly be interesting to verify if these features persist in other models of 
nonabelian vortices, involving different homotopy groups. In the present model, we 
find a Kosterlitz-Thouless type phase structure in the underlying system, i.e. both 
i and j type vortices undergo K-T like phase transitions at some particular critical 
temperatures. What is interesting, and a direct consequence of the nonabelian nature 
of the system, is the connection between these critical temperatures, as evident in 
eq.fTgg). 



23 



Acknowledgements 

The work of S.T. was supported by the Royal Society of Great Britain. 



Appendix A 

Self Energy Modifications 

In this section we wish to show that the term Vm- 

V m = [ d 2 zsm 2 [A(z - z 2 )]d z A(z - z^A{z - z x ) (A.l) 



(see eqs.(4.17) and (4.18) ), can indeed be interpreted as modifications to the divergent 
self energy of test charges placed in the system, as claimed in section 4. In eqn. 
(A.l) Zi, z 2 are the positions of the % type test charge and j type background charge 
respectively. Rather than compute Vm by integrating, we instead consider V 2 Vm since 
this quantity appears in the corresponding P-B equation. One finds 

d Zl d Sl V M = d Zl I d 2 z sin 2 [A(z- z 2 )} 

x d Zl d z A(z - zi)d s A(z - z x ) + d z A(z - z x )d- zx d z A(z - z x ) 
= d Zl fd 2 zsm 2 [A(z-z 2 )][d Zl 5 2 (z-z 1 )-^— + d Zl 5 2 (z-z l )^— 

J ' Z — Z\ z — Zx 

+5 2 (z - Zl )5 2 {z - Zl ) + 7-^—rz] (A.2) 

\Z — Z\\ 

where we remind the reader that A(z — z\) = Im(ln(^ — z\j) 

One can see that all terms in eq.(A.2) are manifestly singular except the last term, 
and hence may be subtracted from the Hamiltonian of the system as discussed in 
section 4. Let us concentrate then on the remaining term, 

sin 2 L4(,-, 2 )] 

\z-z x \ A v ; 
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By a change of variables 

sin 2 [A(z' + LUi)} 



we obtain 



d 2 z 



\z — Z\ 



1 1 

2 " 4 



1 1 
J i 2 ~ 4 



z' +UJj \ f z' + UJj 
Z' + U)i) \Z' + UOi 



z- z 2 \ f z-z 2 



(A.4) 



z - z 2 



,z - z 2 



d 2 z 



\Z — Z\ 



_ I_I 

4 I 2 4 



2Re[zz] + 2Re[z 2 z 2 ] - 4Re[^ 2 ] 



(A.5) 



\Z-Z2\ 

Next, we wish to consider the behaviour of the integral as z — > Z\ . In this respect we 
define 



z 

A 



Z — Z\ 
Z\ - z 2 



r'e 10 



d cos /3 + id sin /3 



(A.6) 



so that we get for the integral in eq.(A.3) 



■2„l 



d 2 r 



„/i 



2 



1 

4 



2 r' 2 cos 26 + 2d 2 cos 2/3-4 r'd cos(6 + (3) 



(A.7) 



(r' + rf) 4 

which has divergent contributions when z —* Z\. In fact the 1/r' term in this expression 
is exactly what we would expect to obtain if we act on the (usual) logarithmically 
divergent self energy of a single vortex in 2-dimensions, with V 2 . 

Having shown that d zl d Sl VM gives rise to divergent terms we can conclude that 
contributions like Vm can be interpreted as self energy modifications to the test charges, 
and hence should be subtracted from the effective hamiltonian of our system. 
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